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ABSTRACT 



We present a detailed and statistical analysis of the column densities and Doppler b parameters of Mgn 
absorbing clouds at redshifts 0.4 < z < 1.2. We draw upon the HIRES/Keck data (Aw ~ 6.6 km s _1 ) 
and Voigt profile (VP) fitting results presented by Churchill & Vogt (Paper I). The sample is comprised 
of 175 clouds from 23 systems along 18 quasar lines of sight. In order to better understand whether 
the inferred physical conditions in the absorbing clouds could be "false" conditions, which can arise due 
to the non-uniqueness inherent in parameterizing complex absorption profiles, we performed extensive 
simulations of the VP analyses presented in this paper. In brief, we find: (1) The Fen and Mgn column 
densities are correlated at the 9 a level. There is a 5 a anti-correlation between the Mgi/Mgn column 
density ratio and the Mgn column density. (2) Power-law fits to the column density distributions for 
Mgn, Fell, and Mgi yielded power-law slopes of ~ —1.6, —1.7, and —2.0, respectively. (3) The observed 
peaks of the Doppler parameter distributions were ~ 5 km s _1 for Mgn and Fen and ~ 7 km s _1 for 
Mgi. The clouds are consistent with being thermally broadened, with temperatures in the 30-40,000 K 
range. (4) A two-component Gaussian model to the velocity two-point correlation function (TPCF) 
yielded velocity dispersions of 54 km s _1 and 166 km s — 1 . The narrow component has roughly twice 
the amplitude of the broader component. The width and amplitude of the broader component decreases 
as equivalent width increases. (5) From photoionization models we find that the column density ratios 
are most consistent with being photoionized by the ultraviolet extragalactic ionizing background, as 
opposed to stellar radiation. Based upon the Mgi to Mgn column density ratios, it appears that at least 
two-phase ionization models are required to explain the data. 

Subject headings: (galaxies:) — quasar absorption lines; galaxies: — halos; galaxies: — kinematics and 
dynamics; galaxies: — ISM 



1. INTRODUCTION 

In Churchill & Vogt (2001, hereafter Paper I), we pre- 
sented 23 Mgn absorption systems observed at high- 
resolution in high signal-to-noise quasar spectra. The 
analysis presented in Paper I was primarily based upon 
direct measurements of the data. In this paper we concen- 
trate on measurements obtained using Voigt profile fitting 
of the data. 

Voigt profile analysis provides a well defined, if simplis- 
tic, parameterization of complex absorption lines. From 
the beginning of the last decade (e.g., Carswell et al. 1991), 
Voigt profile (VP) analysis has been a useful technique 
for statistically documenting the physical conditions of 

1 Based in part on observations obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among Caltech, the 
University of California, and NASA. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. 

2 Center for Gravitational Physics and Geometry 
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absorbing gas in quasar spectra. VP analysis yields the 
number of "clouds" in a given absorption complex, their 
Doppler parameters, and their relative line-of-sight veloc- 
ities. Since the atomic physics is incorporated, the line— of— 
sight integrated column densities and the thermal and/or 
turbulent properties of the gas clouds can be measured. 

The reported statistics of both the high redshift Lya 
forest (Hu et al. 1995; Lu et al. 1996) and Civ absorbers 
(Rauch et al. 1996) have been based upon VP decomposi- 
tion of high resolution (R ~ 45, 000) and signal-to-noise 
ratio {S/N ~ 30 — 80) spectra. A primary assumption 
of VP analysis is that complex absorption profiles can 
be decomposed into a product of individual components, 
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each of which are treated as spatially isolated, isothermal 
clouds intercepting the incoming flux at a unique line- 
of-sight velocity. However, numerical simulations of the 
Lya forest and of Civ systems at high redshift undermine 
this physical picture because the absorbing systems merge 
continuously into a smoothly fluctuating background, of- 
ten are comprised of a range of temperatures, and are 
usually broadened by coherent velocity flows (Haehnelt, 
Stienmetz, & Rauch 1996; Dave et al. 1997). This is not 
altogether surprising for higher redshifts, given that Lya 
and CiV absorption likely arise in large, coherent struc- 
tures that are still coupled to the Hubble flow (e.g., Rauch, 
Haehnelt, & Steinmetz 1997). 

Unblended absorption lines arising in lower ionization 
gas, as traced by the Mgn AA2796, 2803 doublet, are of- 
ten relatively narrow and suggestive of quiescent single 
"clouds" (Churchill ct al. 1999; Churchill & Charlton 1999; 
Churchill & Vogt 2001). As such, the low ionization gas 
more directly associated with the components of galaxies 
may be a better match to the assumptions underlying VP 
analysis. However, Mgll systems also show a wide varia- 
tion in kinematic complexity, including significant blend- 
ing of the profiles arising from the discrete clouds (Petit- 
jean & Bergeron 1990; Churchill & Vogt 2001). It is not 
clear how this particular type of blending could systemat- 
ically skew the statistical results yielded by VP analysis. 
In addition, an unavoidable consequence of VP analysis 
is that the results are sensitive to both the resolution and 
signal-to-noise ratio of the spectra (Churchill 1997; Rauch 
et al. 2002). 

Even with the many caveats mentioned above, the po- 
tential for VP analyses to answer many questions regard- 
ing the physical nature of the absorbing gas remains widely 
recognized. What is the average number of clouds per sys- 
tem? What is the distribution of column densities and 
column density ratios? How does the cloud-cloud veloc- 
ity clustering depend upon equivalent width? Are the line 
broadening mechanisms thermal, or dominated by turbu- 
lent motions in the gas? With the aid of state of the art 
photoionization models we can ask: What are the chemi- 
cal and ionization conditions of the studied ionic species? 
What can be inferred about the nature of the ionizing ra- 
diation? If there is a stellar contribution, what limits can 
be placed on the numbers of stars and their distances from 
the gas clouds? 

For Mgll systems, answers to some of these questions 
have been investigated using data with a resolution of 
30 km s _1 (Petitjean & Bergeron 1990; Bergeron et al. 
1994). Since VP analysis is strongly dependent upon spec- 
tral resolution, it is a logical step to re-examine the sta- 
tistical distributions previously inferred for Mgll systems 
with higher resolution data. In this paper we present sta- 
tistical results from a VP analysis of an ensemble of 23 
Mgn-selected absorbers obtained with the HIRES instru- 
ment (Vogt et al. 1994) on the Keck I telescope. The 
resolution of these data is ~ 6 km s _1 . The systems and 
the VP results were previously presented in Paper I. 

In § 2, we briefly review the data. In § 3, we describe the 
VP analysis of the data, our VP decomposition methodol- 
ogy, and simulations of the analysis for purposes of diag- 
nosing whether or not a given statistical result is induced, 
in part or in full, by the VP technique itself. In § 4, we 



present the distribution of VP column densities, Dopplcr 
parameters, and velocities. We also provide some analysis 
of these distributions. Photoionization models of the VP 
components, or clouds, are discussed in § 5. We include 
an exploration of the possible relative contribution of stel- 
lar flux to that of a standard ultraviolet background. We 
outline the main results of this paper in § 6. 

2. THE DATA 

The acquisition, calibration, and analysis of the data 
has been described in Churchill (1997) and in Paper I. 
However, we provide a brief description here. We em- 
phasize that all systems are selected by the presence of 
Mgll AA2796, 2803 absorption. 

The quasar spectra were obtained with the HIRES in- 
strument (Vogt et al. 1994) on the Keck I telescope. The 
spectral resolution is R — 45, 000. The signal-to-noise ra- 
tio in the continuum near the Mgll absorption profiles is 
typically <~ 30. A total of 23 systems comprise the sam- 
ple of absorbers. The sample includes only those systems 
with Mgll A2796 rest-frame equivalent width, 1^.(2796), 
greater than 0.3 A. The ensemble of systems is considered 
to be representative of a "fair", or unbiased sample, in 
that the distribution of Mgll equivalent widths is consis- 
tent with that of unbiased surveys (e.g., Lanzetta, Turn- 
shek, & Wolfe 1987; Steidcl & Sargent 1992; Churchill et 
al. 1999). 

The sample has been divided into four subsamples by 
their equivalent widths. Sample B has W r = (0.3,0.6] A, 
sample C has W r — (0.6, 1.0] A, sample D has W r > 0.6 A, 
and sample E has W r > 1.0 A. We have used the notation 
"(" to denote inclusive and "]" to denote exclusive. These 
samples are historically motivated (see Stcidel & Sargent 
1992). In Table 1 we present the systems by quasar, ab- 
sorber redshift, and sample membership. Also listed is a 
system identification, or ID. 

The Mgll A2796 profiles are presented in Figure 1. The 
spectra are normalized and converted to rest-frame veloc- 
ity. The systems are grouped by sample membership and 
are ordered by increasing W r . The system ID and equiva- 
lent width are given in each panel. Each system has been 
modeled using Voigt profile (VP) analysis, as described in 
Paper I. The solid curves through the data are the VP 
model spectra and the ticks above the spectra are the VP 
centroids. 

3. VOIGT PROFILE ANALYSIS 

We used the Voigt profile fitting program MINFIT 
(Churchill 1997) to obtain a least-squares fit (LSF), y 2 
minimized model of each absorption line profile. MIN- 
FIT sets up M nonlinear functions in N parameters and 
drives the NETLIB slatec LSF routine DNLS1 (More 1978; 
Hiebert 1980). The N parameters are the redshifts, column 
densities, and Dopplcr b parameters for each VP compo- 
nent. The MINFIT program statistically minimizes the 
number of VP components used to model a given system. 
Numerical convergence is set to the machine precision of 
a Sun Ultra 1. Further details are provided in Churchill 
(1997). 

As discussed in Paper I, we concentrate our analysis on 
the Mgll AA2796, 2803 doublet, the Mgi A2852 transition, 
and the Fen A2344 2374, 2582, 2587, and 2600 multiplet. 
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The atomic data used for the VP analysis are given in Ta- 
ble 2 of Paper I. The column densities, b parameters, and 
cloud velocities output by MINFIT are given in Table 7 
of Paper I and illustrated for all transitions in Figure 2 of 
Paper I. 

The resulting VP model spectra for the Mgll A2796 
transitions are shown in Figure 1 as solid curves through 
the data. Ticks above the spectra give the number of VP 
components and their rest-frame velocities. In Figure 2, 
we show a detail of the fit for both members of the Mgll 
doublet for system S15. For this system, seven individ- 
ual VP components (dotted curves) successfully model the 
complex profile shape. For the illustrated velocity region, 
the reduced chi-square was xl = 1-3 for 418 degrees of 
freedom (including one Mgi and five Fell transitions each 
with 50-60 pixels). 

Since VP analysis yields a non-unique interpretation of 
complex absorption lines, it is important to outline key 
specifics of the fitting procedure and calibration of the re- 
sults if the analysis is to be considered objective. In the 
following, we briefly outline the fitting methodology ap- 
plied in this work, with focus on specifics of the operations 
performed by the MINFIT code. We also describe simula- 
tions of VP analysis that we use to calibrate our statistical 
results. 

3.1. Profile Fitting Methodology 
3.1.1. The Number of VP Components 

MINFIT employs significance testing in order to return 
the minimum number of VP components. Two user- 
defined, run-time parameters, f err and f c , provide the 
maximum fractional error allowed in a VP component and 
the confidence level at which a VP component is consid- 
ered statistically significant, respectively. These parame- 
ters are adjustable, which adds a subjective element to an 
otherwise objective, well defined algorithm. Typical values 
for f err and f c are 1.5 and 0.97, respectively. 

Once a LSF solution is obtained, MINFIT computes er- 
rors in the standard fashion using the diagonal elements of 
the covariance matrix, which are obtained from the inver- 
sion of the curvature matrix. The elements of the curva- 
ture matrix are computed as described in Churchill (1997). 

The fractional error of each VP component is the 
quadrature sum of the fractional errors in the three VP 
parameters, N, b, and z. When the fractional error of a 
component exceeds f err , the component is dropped. The 
LSF and error computations are iterated until no compo- 
nents have fractional errors exceeding f err . 

At this stage the significance level is determined for the 
VP component with the largest fractional error. The tar- 
get VP component is dropped and a new LSF solution is 
obtained. Then, an F-distribution test is performed be- 
tween the two LSF solutions. If the VP component was 
significant with a confidence level of f c or better, than the 
original LSF model is restored. If not, the statistical sig- 
nificance is determined for the VP component with the 
largest fractional error in the new model. This process 
is repeated until all VP components are significant at a 
minimum confidence level of f c . 



3.1.2. Redshifts 

It was assumed that the redshifts, z, and therefore rest- 
frame velocities, of each VP component are aligned for 
all ions. VP components were defined by the presence of 
Mgll absorption features even if there were no Fen or Mgi 
absorption features at the same velocity. 

We note that the ionization potential of Mgi is below 
the hydrogen edge. The adopted methodology of enforc- 
ing an aligned velocity for all ions may not be sound if the 
gas has an ionization structure that correlates with line of 
sight velocity. That is, Mgi could, in principle, arise from a 
physically and/or kincmatically distinct part of the clouds 
in comparison to Mgll and Fell. 

3.1.3. Column Densities 

For each VP component, the column densities, N, of 
different ions are determined independent of one another. 
When multiple transitions of a given ion are available, then 
the increased number of constraints yield a more robust 
measurement. For Mgi there is only the singlet. For Mgll 
there is always the resonance doublet and for Fen there 
are often three to five transitions of the multiplet available. 
As such, Mgi is the least constrained ion. 

If neither Mgi nor Fen had detectable absorption at the 
same velocity as a Mgll component, then an upper limit 
on the Mgi or Fen column density is quoted for the com- 
ponent. This limit is taken from the 3 a equivalent width 
threshold for an unresolved absorption line at the expected 
location of the component (see Churchill et al. 2000a). 

3.1.4. Doppler Parameters 

The line broadening, or b parameter, of a VP component 
can be written as the convolution of a thermal component 
and a turbulent component. In order to preserve the Voigt 
profile formalism, one is required, perhaps quite unrealis- 
tically, to parameterize the turbulence with a Gaussian 
distribution, yielding b 2 = b 2 therm + b 2 turb . Whereas b t h erm 
scales with the inverse square root of the ion mass, bt U rb is 
independent of ion mass and will have the same magnitude 
for all ions. 

There are three possible treatments of the b parame- 
ters. In each VP component, one can: (1) measure the b 
parameter independently for each ion with no assumptions 
about thermal scaling or level of turbulence, (2) force the 
b parameter to be equal for all ions, which holds if the tur- 
bulent component dominated, or (3) force the turbulent 
component to be some fraction of the measured b parame- 
ter (zero for pure thermal). We assume the first approach, 
so that we can fully explore the line broadening mecha- 
nism. 

3.2. Calibrating Our VP Analysis 

The only way to calibrate VP analysis is to quantify the 
degree to which the fitted parameters reflect the "true" 
cloud properties in cases where they are known a priori. 
This requires simulations for which the VP parameters ob- 
tained from analysis of synthetic spectra can be directly 
compared to the VP parameters used to generate these 
spectra. 

We have run a series of comprehensive simulations of 
Mgll absorbers for varying degrees of line blending (i.e. 



4 



CHURCHILL, VOGT, & CHARLTON. 



complexity of absorption profiles) and for a range of spec- 
tral signal-to noise ratios (i.e. equivalent width detection 
threshold) commensurate with the data. These simula- 
tions have been described in detail in Churchill (1997), 
and serve to quantify any possible bias, systematic, or 
random errors in the VP results to the data. For each 
analysis using VP fits to the observational data, we per- 
formed an identical analysis using VP fits to the synthetic 
spectra. In this way, the systematics or significance of any 
observational statistical result can be quantified and any 
systematic skew in the functional form of a distribution 
can be estimated. 

The simulations required input column density, 6 pa- 
rameter, and velocity clustering distributions. We used 
the data to guide our "best" choices of these distribution 
functions by comparing (1) the distribution of pixel flux 
decrements, (2) the distribution of equivalent widths of 
the kinematic subsystems 3 , (3) the distribution of velocity 
widths of the kinematic subsystems, and (4) the observed 
VP component column density, 6 parameter, and velocity 
distributions. 

The adopted distributions for the simulations were ob- 
tained by constructing a set of 500 simulated absorp- 
tion systems and performing VP analysis on these spec- 
tra. We then performed statistical tests comparing the 
simulation distributions and the observed distributions 
to eliminate those which were not statistically consistent 
with data (also see Charlton & Churchill 1998). For the 
Mgn column densities, we chose a power law distribution, 
f(N) oc N~ s with S = 1.6. The b parameters distribu- 
tion was drawn from a Gaussian with centroid 4 km s _1 , a 
width of o~ == 1.5 km s _1 and a lower cut-off of 1.5 km s _1 . 
The lower cut off was required to match the small b por- 
tion of the observed distribution. For the VP component 
velocity clustering, we chose a Gaussian distribution with 
(j = 45 km s _1 . 

We also simulated the Fen multiplet. We enforced ther- 
mal broadening for all VP components, i.e. the b parame- 
ters scaled as 6(Feii) = {m(Mg)/m(Fe)} 1/2 6(Mgii), where 
m is the atomic mass. We also assumed a linear relation- 
ship between log TV (Fen) and logTV(Mgn) for all compo- 
nents. This relation, log N (Fen) = logiV(Mgll) - 0.3, is 
based upon preliminary VP fitting to the sample presented 
in Churchill (1997). 

Similar simulations for VP analysis calibration have 
been undertaken by Lu et al. (1996) and by Hu et al. (1995) 
for the Lya forest. The significant difference between the 
simulation performed here and those for the Lya forest are 
that (1) the hydrogen absorption-selected clouds were as- 
sumed to be distributed cosmologically, whereas the Mgn 
absorbing clouds are clustered in kinematic systems, and 
(2) the VP analysis of the Mgn systems included simulta- 
neous fitting to the Mgn doublet and Fen multiplet. 

4. RESULTS OF VOIGT PROFILE ANALYSIS 

4.1. Number of VP Components (i.e. Clouds) 

In Figure 3a we show the distribution of the number of 
VP components (i.e. clouds). The dotted-line histogram 
represents the full sample of 23 systems. The shaded his- 
togram represents a subset of the full sample in which 

3 Kinematic subsystems were defined in Paper I. They are the individual absorption features within ±500 km s _1 of the Mgn absorption 
centroid that are separated by a minimum of three pixels (a resolution element) of continuum flux. 



the very large iV(Hi) systems (S2, S5, and S7) have been 
removed. These damped Lya absorbers (DLAs) and "Hi- 
rich" systems (Churchill et al. 2000b) have been removed 
in light of their characteristic black-bottomed Mgn ab- 
sorption profile morphologies; it is difficult to constrain 
the number of clouds in these systems (they are system- 
atically underestimated unless there are good constraints 
from weak Fell transitions). For ease of discussion, we will 
refer to this subset of our 23 systems as the "non-DLA" 
subset. 

For both the full sample and the non-DLA subset, there 
is an average of 7.7 clouds per absorber. The mode of 
the binned distribution is eight clouds. In the simulations 
of complex, blended profiles, ~ 30% of the inputted VP 
components were not recovered. These simulations also 
revealed that the percentage of unrecovered components 
decreases with decreasing signal-to-noise ratio, but that 
the width of the recovered distribution becomes larger as 
the signal-to-noise ratio increases. To the extent that VP 
components (such as those in our models) give rise to the 
observed absorption profiles, the distribution plotted in 
Figure 3a would require a correction such that the aver- 
age and mode is elevated to 10-11 clouds per system. 

In Figure 36, we show the rest-frame equivalent width 
of Mgn A2796 versus the number of VP components, N c i. 
The unfilled data points are the DLA and Hi-rich sys- 
tems. As has been reported previously (Petitjean & Berg- 
eron 1990; Churchill et al. 1996; Churchill 1997), there is a 
very strong correlation between the equivalent width and 
the number of VP components. A linear LSF, including 
errors in the equivalent widths, was performed on the data 
presented in Figure 36. The dotted line is the fit to the 
full sample, which yielded a slope of 0.058 ±0.004 A/cloud 
and an equivalent width of 0.28 ± 0.03 A for N c i = 1 (the 
intercept at N c i = has no physical meaning). 

A fit to the non-DLA subset yielded a statistically iden- 
tical slope with a slightly smaller value of 0.23 ±0.03 A for 
N c i — 1. This fit is shown as the solid line on Figure 36. 
Thus, at least in this sample of 23 systems, the three DLA 
and Hi-rich systems serve only to slightly increase the 
intercept to the fit and do not skew the slope. We also 
examined whether the data point at N c i = 18 was dom- 
inating the fit results by removing it and refitting to the 
non-DLA subset of data; the resulting fit was statistically 
identical to the fit to the full sample. 

The slopes of the linear relationship presented here 
are slightly smaller than both the slopes of 0.070 ± 
0.004 A/cloud reported by Churchill et al. (1996) for a 
smaller sample of 15 systems and 0.076 ± 0.004 A/cloud 
reported by Churchill (1997) for a sample of 36 systems. 
This difference is most likely due to the fact that in this 
work we did not enforce the fit to pass through the origin, 
whereas the fit was forced through the origin in the previ- 
ous work. We stress that the slope is strongly dependent 
upon the spectral resolution of the data. For example, Pe- 
titjean & Bergeron (1990) obtained a much steeper slope 
of 0.35 A/cloud for a large sample observed at a resolution 
of ~ 30 km s _1 . 
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4.2. Column Densities 

In Figure 4, we present Fell and Mgl column densities 
vs. Mgn column densities. Only components for which 
the errors are less than 0.7 dex are plotted. Upper limits 
on the Fen and Mgl column densities are shown as open 
points with downward pointing arrows and measured val- 
ues are shown as solid points with 1 a error bars. 

A(Feii) is strongly correlated with A(Mgn). A 
Spearman-Kendall test yields a 9 a significance. As- 
suming a linear relationship of the form log AT (X) = 
a log A(Mgn) + b, we determined the best fit parameters 
using a \ 2 minimization that allowed for errors in the Fen, 
Mgl, and Mgn column densities. The best fit parame- 
ters are a = 0.73 ± 0.06 and b = 3.0 ± 0.8 for Fen and 
a = 0.45 ± 0.05 and b = 5.2 ± 0.6 for Mgl. Limits have 
been included in the fitting. 

The reduced \ 2 is 30 for Fen and is 27 for Mgl. Fits 
to the Fen vs. Mgn and Mgl vs. Mgn column densities 
obtained from the simulations yielded \ 2 — 10- Thus, the 
variance in the data is roughly a factor of three greater 
than the variance due to statistical scatter inherent in the 
VP fitting. 

In Figure 5, we present Fell to Mgn and Mgl to Mgn 
column density ratios vs. Mgn column densities. The data 
presented are the same VP components shown in Fig- 
ures 4a and 46. Upper limits on the ratios are shown 
as open points with downward pointing arrows and mea- 
sured values are shown as solid points. The error bars 
are suppressed, but they are effectively identical to those 
presented in Figure 4. 

The distribution of A(Feii)/A(Mgn) visually appears 
to be a slightly decreasing function of A(Mgn), but is sta- 
tistically consistent with no dependence upon A(Mgn). 
This would indicate that A(Fell)/A(Mgn) is fairly in- 
dependent of or only slightly decreasing with increasing 
A(Mgn). The greater than 1 dex spread in the ratio would 
indicate a range of [a/Fe] abundance patterns, dust deple- 
tion factors, ionization conditions, or a combination of all 
three. 

There is a significant (5 er) anti-correlation between 
A(Mgl)/A(Mgn) and A(Mgn) over a 2.5 order of mag- 
nitude range of the ratio. For the upper limits, taken by 
themselves, this trend is an artifact of the limiting de- 
tectable A (Mgl), which is effectively a constant indepen- 
dent of A(Mgll). However, for the measured ratios, the 
anti-correlation is real and cannot easily be explained as 
an artifact of the VP analysis of the data. A linear LSF 
to the detections only (filled data points of lower panel in 
Figure 5) yielded a slope of a = —0.73 ± 0.06 and inter- 
cept b = 7.6 ± 0.3. The linear fit and its uncertainty are 
superimposed on the data. 

Further discussion of the column density ratios is pre- 
sented in § 5 in the context of photoionization models. 

4.3. Column Density Distribution 

VP analysis of simulated spectra with comparable 
signal-to-noise ratios revealed that the input column den- 
sity distribution is recovered within the 1 er errors when 
an appropriate lower-limit, cut-off column density is ap- 
plied to a maximum likelihood fit. This cut off is the 
column density below which the detection completeness 
drops rapidly. 



The simulations further revealed that the completeness 
level is sensitive to the kinematic clustering of components 
due to line blending. In regions of significant blending 
the 90% completeness levels are log A(Mgn) = 12.4 cm~ 2 
and log A(Feii) = 12.2 cm~ 2 (we did not examine blend- 
ing effects for Mgl). The 90% completeness levels for un- 
blended lines are log A(Mgn) = 11.6 cm~ 2 , log AT (Fell) = 
11.8 cm~ 2 , and logA(Mgi) = 11.4 cm~ 2 . Thus, there is 
a column density range of "partial completeness" , a tran- 
sition from incompleteness due to line blending to incom- 
pleteness due to finite signal-to-noise ratio. 

In Figure 6, the column density distribution functions 
are shown for Mgn, Fen, and Mgl. The data arc binned 
for purpose of presentation. The shaded regions show 
the column density ranges of partial completeness. The 
data appear to follow a power-law distribution, with the 
turnover at small column densities due to incompleteness. 
At HIRES/Keck resolution, an unresolved Mgn line satu- 
rates at log A ~ 13.5 cm~ 2 and this explains the increased 
scatter in the number density of Mgil VP components with 
column densities greater than this value (see Figure 4.3 of 
Churchill 1997). 

For the column density distribution, we assumed a 
power law 

f(N) = CN~ S , (1) 

where /(A) is the number of clouds with column density 
A per unit column density. The total number of clouds, 
m, normalizes /(A) over the observed data according to 

/ f(N)dN = m. (2) 

We used the maximum likelihood method (e.g., Tytler et 
al. 1987; Lanzetta, Turnshek, & Wolfe 1987) to obtain the 
power-law slope and normalization. The method is per- 
formed on the unbinned data. In Figure 6, the solid lines 
through the data are the maximum likelihood solutions, 
which are 

Mgn: S = 1.59 ±0.05 C = (9.85 ± 0.16) x 10 8 , (3) 

Fen: 5 = 1.69 ±0.07 C = (1.01 ± 0.22) x 10 10 , (4) 

Mgl : 5 = 2.02 ±0.28 C = (6.21 ± 0.79) x 10 12 . (5) 

For this analysis, we used the cut-off column densities ap- 
propriate for blended lines, which yielded m — 161, 107, 
and 71 for Mgil, Fen, and Mgl, respectively. 

Using the parameterization given in Equation 1, Petit- 
jean & Bergeron (1990) found a 5 = 1.0±0.1 for Mgn over 
the range 12.0 < logA(Mgn) < 14.3 cm" 2 . The resolu- 
tion of their data was ~ 30 km s _1 , and it is likely that 
their best fit slope is affected by unresolved blending. The 
expected trend is that the number of large column den- 
sity components would decrease, being redistributed into 
a greater number of smaller column density components as 
the resolution is increased. This would result in a steeper 
slope, as we have found in this work. 

The slopes of the Mgn and Fen column density dis- 
tributions are very similar, and roughly consistent with 
one another. However, the slope of Mgl distribution is 
significantly steeper. These findings are consistent with 
the distribution of column density ratios presented in 
Figure 5. There is no correlation of A(Feii)/A(Mgn) 
with A(Mgn), yet there is a significant anti-correlation 
of A(Mgi)/A(Mgn) with A(Mgn). 
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4.4. Doppler Parameter Distribution 

The binned distribution of Mgn, Fen, and Mgi Doppler 
b parameters is presented in Figure 7. The median Doppler 
parameters and the standard deviations of the distribu- 
tions are (6) = 5.4 ± 4.3, 5.1 ± 4.1, and 7.7 ±5.1, for Mgil, 
Fen, and Mgi, respectively. We note that, statistically, 
the relative distributions of Mgil and Mgi are unphysical; 
it is difficult to understand how Mgi could have broader 
lines than Mgil on average. 

Based upon simulations, the observed distribution peak 
for Mgil could be shifted to a larger value by ~ 1-2 km s _1 
compared to the "true" distribution. The true Mgil dis- 
tribution is likely peaked at b = 3-4 km s -1 . Simulations 
further indicate that the true distribution may have zero 
to few clouds with b < 1 km s _1 . A lower cut-off b param- 
eter of 1.5 km s _1 is required to not overproduce the low- 
est bin. However, we note that this could be a resolution 
effect, given that lines are unresolved for b < 2.46 km s" 1 . 

In 30 km s _1 resolution data, Petitjean & Bergeron 
(1990) measured a Mgil b distribution peaking at 10 < b < 
15 km s , with a decaying tail at b > 20 km s^ 1 . They 
suggest that non-thermal or turbulent motions within the 
individual clouds are implied because such large b values 
correspond to temperatures in excess of 10 5 K, a tempera- 
ture too high for Mgil to survive if photoionized. Our re- 
sults are suggestive of clouds that are less turbulent and/or 
cooler, with temperatures between 30-40, 000 K. In our 
data, all VP components with b > 10 km s -1 occur in 
broad, partially or fully saturated profiles. 

Simulations revealed that the majority of the large b pa- 
rameters forming the extended tail to the distribution is 
an artifact of component blending. Approximately 70% of 
the input VP components were recovered; that is, three of 
ten components were typically lost because of blending 4 . 
As such, some components may have large b parameters 
in cases where two components may have been responsi- 
ble for the absorption profile, but could not be formally 
constrained by the fitting algorithm. 



4.5. Thermal and Turbulence Broadening 

The b parameter of a given ion can have both a ther- 
mal and non-thermal component. This latter contribution 
to the line broadening could be due to internal turbulent 
motion or systematic, bulk motion. Assuming the non- 
thermal component can be parameterized as a Gaussian 
distribution preserves the Voigt profile formalism. We note 
that a Gaussian non-thermal motion in an isolated cloud is 
difficult to understand, since it implies a mixing, or quasi- 
convective process that would be expected to decay away 
in a cloud dynamical time. 

The total b parameter of ion i is written, 



h 2 —h 2 ■+■ h 2 

u i — u i,therm < u , 



turbi 



(6) 



where bitherm — 

^2kTjm l is the thermal component, 
which scales with the ion mass m, and b tur b is the non- 
thermal component, which is equal for all ions. The ther- 
mal and non-thermal b parameters for two ions of different 
mass are determined by solving two equation with two un- 

4 This result is sensitive to the signal-to-noise ratio of the spectra, 
data. 



knowns, which gives 
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where bi and bj are the measured b parameters for the 
lighter and heavier ions, respectively. This technique for 
exploring the line broadening mechanism has been applied 
using Civ and Siiv at z ~ 2.5 by Rauch et al. (1996). 

In Figure 8, we present the results of deconvolving the 
thermal and non-thermal components to the line broaden- 
ing of Mgil and Fen components. The selected data have 
b parameters with fractional errors less than 50% for both 
Mgil and Fell. The thermal component is plotted along 
the vertical and the total b parameter is plotted along the 
horizontal axis. The solid curves are lines of constant non- 
thermal motion. The temperature scale for the thermal 
component is given on the right hand vertical axis. 

Clouds that lie above the bturb = km s _1 curve have 
unphysical Mgil to Fell b parameter ratios (recall that the 
b parameters were fit assuming no constraining relation- 
ship between different ions). The mean turbulent compo- 
nent of the subsample is 0.7 ± 0.3 km s _1 . VP analysis of 
simulated spectra in which the clouds are assumed to be 
thermally broadened, i.e. have b tU rb = km s _1 , yielded a 
mean b tur b of ~ 1 km s _1 . This value is consistent with the 
value obtained for the data and is suggestive of thermally 
broadened clouds. 

Furthermore, the results from the simulated spectra 
showed the same level of scatter as the observed data plot- 
ted on Figure 8. As such, it is difficult to claim a non- 
thermal component in individual systems on a case-by- 
case basis. Based upon the simulations, it is found that 
the scatter decreases with increasing signal-to-noise ratio. 
We also note that even though the multiple Fen transi- 
tions offer greater constraints on the Fen b parameter, the 
weaker strength Fen data often have signal-to-noise ratios 
a factor of 2-3 lower than the Mgil data. 

The simulation results and the spread of the data about 
the bturb = km s _1 curve provide no compelling statis- 
tical evidence for a non-thermal component to the line 
broadening of Mgil absorbing gas. From thermal compo- 
nents, the cloud temperatures are inferred to be in the 
range 30, 000-150, 000 K. Based upon the photoionization 
modeling, it is difficult to understand temperatures in ex- 
cess of 60, 000 K. Again, we note that the clouds with the 
largest b parameters are due to the blending of compo- 
nents. 

4.6. Cloud-Cloud Velocity Clustering 

We computed the velocity two-point clustering func- 
tion (TPCF) for the VP components obtained with our 
HIRES/Keck spectra (these data are listed in Table 7 of 
Paper I). As computed here, the TPCF, P(Av), is the 
probability that any randomly selected pair of VP compo- 
nents in a system will have a velocity splitting Av. 

The quoted results are for signal-to-noise ratios commensurate with the 
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Using low resolution spectra, Sargent, Steidel & Bok- 
senberg (1988) used the TPCF to show that Mgn sys- 
tems clustered on velocity scales of Av < 200 km s _1 . 
This was taken as evidence that Mgn systems cluster 
like field galaxies. Using moderate resolution spectra 
(~ 30 km s _1 ), Petitjean & Bergeron (1990) found that 
a two component Gaussian model with a\ = 80 and 
er 2 = 390 km s- 1 provided a good fit to the Mgn TPCF. 
They interpreted the smaller width to reflect the kinemat- 
ics of clouds bound within galactic halos and the broader 
width to reflect the kinematics of galaxy pairs in the field. 

The TPCF for the full sample is presented in the up- 
per panel of Figure 9. The TPCF increases sharply 
for Av < 100 km s _1 and has an extended tail out to 
~ 400 km s _1 . Following Petitjean & Bergeron (1990), we 
parameterized the TPCF as a two-component Gaussian, 

P(A,)=A 1 exp^j+A 2 exp^j ) (9) 

where the A and a are the amplitudes and Gaussian 
widths, respectively. We obtained a narrow width of 
ui = 56 km s~ x and a broader width of a 2 = 166 km s _1 
for the full sample. The relative amplitudes differ by a 
factor of two, with the narrow component dominating. 

We examined whether the extended tail of the distri- 
bution is dominated by the three kincmatically complex 
systems, SI, S9, and S12 (see Figure 1). Not only do these 
systems have a large kinematic spread, they have roughly 
twice the Lya absorption strengths and very strong Civ 
absorption compared to "classic" Mgn systems; they are 
classified as "double systems" (Churchill et al. 2000b). 

If the double systems are removed from the sample and 
the TPCF recomputed, we obtain the results presented in 
the lower panel of Figure 9. For this "No Doubles" sam- 
ple, we obtained a narrow width of o\ = 60 km s _1 and a 
broader width of o~2 = 151 km s _1 . The component widths 
are not substantially different from the full sample TPCF. 
However, the amplitude of the broad component is a factor 
of five smaller than that of the narrow component. This 
indicates that the double systems strongly govern the am- 
plitude of the broad component to the full-sample TPCF, 
but that the component still results because of "moderate- 
" and "high velocity" kinematic subsystems 5 in the "clas- 
sic" Mgn absorbers (Churchill et al. 2000b). 

4.7. Evolutionary Clues from Kinematics 

It is difficult to apply a straight-forward, physical inter- 
pretation to each of the two components of the TPCF. 
However, based upon absorber models (Churchill 1997; 
Charlton & Churchill 1998), the narrow component width 
may be proportional to the vertical velocity dispersion of 
gas in galaxies disks (averaged over random inclination 
angles). This is also consistent with the notion that the 
dominant kinematic subsystems (see Paper I), could arise 
in systematically rotating disks (e.g., Lanzetta & Bowen 
1992; Charlton & Churchill 1998; Churchill & Vogt 2001). 

The data directly constraining this scenario are few in 
number; in four out of five edge-on spiral galaxies host- 
ing TU r (2796) > 0.3 A absorption, the gas kinematics di- 
rectly trace the outer extensions of the stellar disk rotation 
curves (Steidel et al. 2002). However, the data arc difficult 

5 Kinematic subsystems are defined in Paper I. 



to understand in terms of a simple disk model. Counter 
examples are also few; Lanzetta et al. (1997) and Ellison, 
Mallen-Ornelas, & Sawicki (2002) have each reported an 
example inconsistent with a rotating disk scenario. 

Clues to the nature of the broader component to the 
TPCF might be connected to the rcdshift evolution of 
W r (2796). Steidel & Sargent (1992) showed that as the 
minimum equivalent width of a Mgn-absorber sample is 
increased, the redshift number density evolution becomes 
stronger. That is, the strongest systems evolve away most 
rapidly with decreasing redshift (increasing cosmic time) 
from z ~ 2 to z ~ 0.3. The underlying physical pro- 
cess governing this evolution remains a mystery, though it 
has been attributed to either a reduction in the mean col- 
umn density of the clouds with time (chemical, ionization, 
and/or gas mass evolution) or a reduction in the velocity 
spread of the clouds with time (kinematic evolution), or 
both. 

To explore trends in the kinematics with equivalent 
width, we computed a separate TPCF for samples B, C, D, 
and E, with and without double systems. The results are 
listed in Table 2. From these distributions, we see that the 
nature of the cloud-cloud clustering is strongly connected 
to the equivalent width. 

Including double systems, the velocity dispersion of the 
narrow component, <7i, does not change significantly with 
increasing W r (2796) , whereas the broad component veloc- 
ity dispersion, 02, decreases and its relative amplitude, 
A2/A1, increases. If our sample is taken to be represen- 
tative of the population of Mgn absorbers, this trend in 
the TPCF would suggest that the frequency of highest 
velocity clouds decreases as equivalent width is increased 
(note that the exclusion of the double systems makes no 
change to the results for sample B.) In other words, in 
the largest equivalent width systems, the cloud velocities 
are distributed more uniformly throughout the full veloc- 
ity interval of the profile and clouds at extreme velocities 
are less frequent. 

Excluding double systems, the extended tail in the 
TPCF vanishes for samples D and E. Apparently there are 
two types of kinematics characterizing the very strongest 
equivalent width systems, and they are extreme to one an- 
other. The first type is exemplified by systems SI and S9 
(double systems), where the profile morphology is charac- 
terized by a large velocity spread of many clouds repre- 
senting a wide range of optical depths. Systems S5 and 
S7 (which are DLAs or near-DLAs) are examples of the 
second type. The Mgn profile morphology of DLAs and 
near-DLAs, in the studied redshift regime, are character- 
ized by fully saturated, black-bottom absorption with a 
moderate kinematic spread of ~ 60 km s _1 (Churchill et 
al. 2000b). 

These findings may provide insight into the nature of 
the differential equivalent width evolution quantified by 
Steidel & Sargent (1992). Since the redshift number den- 
sity of DLAs apparently does not evolve over the redshift 
range 0.5 to 2.0 (Rao & Turnshek 2000), those large equiv- 
alent width Mgn systems arising in DLAs would also not 
evolve over this redshift interval. Therefore, if the kine- 
matics of the largest equivalent width systems do fall into 
the above two categories, and one type arises primarily in 
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DLAs, then it must be the double systems undergoing the 
most rapid evolution. 

This would suggest that Mgn absorption cross section 
of the kinematically higher velocity material is decreasing 
with time. However, moderate- and high-velocity kine- 
matic subsystems with weak absorption do not contribute 
significantly to the total equivalent width in a system. 
Thus, the evolution cannot be dominated by a decrease 
in the frequency of weak subsystems in double systems 
like SI and S12. It must be an evolution in systems more 
like S9, where the optical depth is highly variable and rela- 
tively large across the entire kinematic range of absorption. 
At even larger equivalent widths, ~ 2 A and greater, the 
evolution is likely dominated by systems with kinematics 
similar to those studied by Bond et al. (2001). A sample 
of high resolution Mgn systems at z > 1 will be required 
to ultimately address these issues. 

4.8. Column Densities and Kinematics 

In Figures lOa-lOc, we present the logarithmic column 
densities for Mgn, Fen, and Mgl. Solid points are mea- 
surements and open points are upper limits. The anti- 
correlations with velocity are induced by the definitions of 
the velocity zero point. The zero points are defined as the 
optical depth median of the absorption profiles and virtu- 
ally all of them are aligned with the centroid of a dominant 
kinematic subsystem (see Paper I). The dispersion of the 
column densities decreases with velocity, though it is not 
clear whether the very small dispersion for v > 200 km s _1 
is due to small number statistics. 

For logiV(Mgii) < 12.5 cm~ 2 , the distribution of Mgn 
column densities is quite flat with velocity, which indicates 
that low column density clouds are present at all veloci- 
ties. For these clouds, Fen and Mgl are often below the 
detection threshold of the data. This explains the pres- 
ence of upper limits for Fell and Mgl at all velocities, and 
most obviously at high velocities. 

For our sample, we can at least say that it is rare to 
find large column density clouds at large velocities from 
the dominant kinematic subsystem (also see Paper I). As 
trivial as this observation may seem, it has implications for 
the systematics in the absorption kinematics and therefore 
possibly the absorber geometry. Apparently, finding two 
large kinematic subsystems in a single absorber is rare. 
These subsystems are comprised of several blended com- 
ponents, which is suggestive of gas that is part of a coher- 
ent structure with low velocity dispersion, i.e. systematic 
kinematics. Since there is only a single dominant kine- 
matic subsystem per absorber and if its kinematics are 
systematic, the geometry is likely to be quasi-planar with 
velocities primarily parallel to the plane. 

In Figures lOd and lOe, we present log iV(Feii)/iV(Mgii) 
and log V(Mgi)/V(Mgn) vs. cloud velocity. Interest- 
ingly, the dispersion of logiV(Feii)/iV(Mgii) with veloc- 
ity is fairly constant with a range from —1 to +0.2 out 
to v ~ 200 km s _1 . This would suggest that the range 
of physical conditions (abundance pattern, dust depletion, 
ionization conditions) governing this ratio does not change 
dramatically with kinematics. If, as suggested above, 
the main kinematic subsystems are single coherent struc- 
tures, then it may be that the higher velocity clouds are 
"cloudlettes" surrounding this structure. If so, then these 



"cloudlettes" apparently do not differ greatly in their Fen 
to Mgn conditions. Beyond v ~ 200 km s _1 we can only 
comment that there are no clouds in our sample with ele- 
vated V(Fell)/V(Mgn). 

Since Mgl absorption is present primarily at velocities 
less than ~ 100 km s _1 , it is less clear if the dispersion of 
log A r (Mgi)/A r (Mgn) with velocity is decreasing at larger 
velocities. In fact, 65% of the detected Mgl arises within 
v < 30 km s -1 and 80% arises with v < 60 km s _1 . 
The range of W(Mgl)/W(Mgll) for v < 60 km s"\ how- 
ever, is four orders of magnitude (excluding the point at 
log V(Mgi)/V(Mgn) = 0). This suggests that the con- 
ditions where Mgl arises are probably quite varied. This 
is likely an ionization condition effect. It could be due 
to differing cloud densities, ionizing flux, or to multiphase 
ionization conditions (see next section). 

5. PHOTOIONIZATION MODELS 

A grid of Cloudy 90 (Ferland 1998) photoionization 
models were constructed in order to place constraints on 
the chemical and ionization conditions of the VP compo- 
nents. Several researchers have applied photoionization 
models to Mgll-selected absorption systems (e.g., Berg- 
eron & Stasinska 1986; Bechtold, Green, & York 1987; 
Bergeron et al. 1994; Churchill & Charlton 1999). A use- 
ful exposition on the methods used here has been written 
by Steidel (1990). 

In brief, the clouds are assumed to be plane-parallel, 
constant density "slabs" with an external ionizing flux in- 
cident upon one face of the cloud. For this configuration, 
the ionization condition is quantified using the ionization 
parameter, U = n^/na, where n 7 is the number density of 
hydrogen ionizing photons and nn is the number density 
of hydrogen atoms. The quantity n 7 is dependent upon 
both the intensity (amplitude at 1 Ryd) and shape of the 
ionizing spectrum. The clouds are given a pre-defined 
neutral hydrogen column density, 7V(Hi). This results in 
well-defined clouds that are indexed by their U and iV(Hi) 
values. We produced a grid of model clouds over the range 
-5 < \ogU < -1 and 15 < logJV(Hl) < 20 cm" 3 at half 
dex intervals in both quantities. 

In these photoionization models, the column density ra- 
tios of various ions are effectively independent of metal- 
licity. Only when metallicity approaches the solar value, 
resulting in increased cooling rates, are the ratios modi- 
fied from those in lower metallicity clouds. We use l/10th 
solar metallicity [Z = —1) for our models. We also use a 
solar abundance pattern. 

5.1. Ionizing Spectrum 

The spectral energy distribution (SED) and intensity 
of the ionizing flux is most uncertain. We explored both 
an extragalactic ultraviolet background (UVB) source and 
"local" stellar sources. For the UVB, we used the spec- 
tra of Haardt & Madau (1996). These spectra have only 
slightly different SEDs for z = 1 and z = 0.5 with normal- 
izations at 1 Ryd of logf/„ = —5.2 and —5.6 erg s' 1 cm -1 , 
respectively. We find that there is negligible difference in 
the column density ratios of Fen and Mgl to Mgil for 
cloud models ionized by these two UVB spectra. For the 
following discussions, we use the z = 1 SED. 

Assuming no shielding of the UVB, stellar radiation 



VP ANALYSIS OF Mgn ABSORBERS 



9 



must dominate over the UVB before the stellar SEDs be- 
gin to alter column density ratios in cloud models. The 
ionization potentials to destroy Mgn and Fell are slightly 
greater than 1 Ryd, whereas the ionization potential to 
destroy Mgi is slightly below. As shown in Appendix B 
of Churchill & Le Brun (1998), only O and B stars are 
capable of modifying the UVB in this energy range for 
astrophysically reasonable stellar number densities. They 
also reported that their results were insensitive to stel- 
lar surface gravity. Thus, we limited our exploration to 
stars with T = 30, 000 K, to represent B stars, and with 
T = 50, 000 K, to represent O stars. Following Churchill 
& Le Brun (1998), we used the solar metallicity stellar 
models produced by Kuruzc (1991). 

As with the UVB grids, we produced a grid of model 
clouds over the range —5 < \ogU < —1 and 15 < 
log AT (Hi) < 20 cm~ 3 at half dex intervals in both quan- 
tities. We also varied the log vf v of stars at 1 Ryd over 
the range —7 to —3 erg s _1 cm~ 2 at 1 dex intervals. We 
produced a separate grid for B and O stars, respectively. 
For the stellar SEDs, there is a sharp edge at 1 Ryd with a 
relatively flat energy distribution at energies just greater 
than 1 Ryd. Thus, for subsequent discussion, we focus on 
the normalization of the stellar grids at 1.2 Ryd, which 
is very close to the ionization potentials to destroy Mgn 
and Fen. For reference, the z = 1 UVB at 1.2 Ryd has 
logvf v = —5.3 erg s _1 cm~ 2 . 

5.2. Limits on Stars and Their Distances 

Following the formalism of Churchill & Le Brun (1998), 
we derive the minimum number of O, B I, and B V stars 
in a confined region at a distance of R kiloparsecs, that 
are required to influence the ionization conditions to be 

IogJV.(O)=0.1 + 2IogiJ, (10) 

log A* (B I) = 0.0 + 2 log R, (11) 

log AT* (B V) = 0.8 + 2 log i?, (12) 

at z = 1. The constants arc reduced by 0.4 for z = 0.5. 
These constants were computed using the flux level of the 
UVB SED at 1.2 Ryd, effectively at the ionization po- 
tentials of Mgn and Fen. To recompute A* for a known 
stellar flux, uf v at 1.2 Ryd, one simply adds 5.3 + \ogvf v . 

We have assumed that the photon escape fraction from 
the stellar environment is 100%. Under the assumption 
of a grey opacity, log A* would be increased by | log f esc | , 
where < f esc < 1. In Figure 11, we plot log A* vs. logi? 
for the three stellar types at z = 1. We apply these results 
below. 

5.3. Fen to Mgn Ratio 

In Figure 12, we show results of the photoionization 
models for the Fen to Mgn column density ratios over- 
plotted on the data. Six panels are shown. The upper 
left panel shows the SEDs for the assumed ionizing con- 
tinuum. The thick solid curve is the Haardt & Madau 
UVB at z = 1. The thinner solid curves are the SEDs of 
T = 30, 000 (B-type) stars for flux levels \og{vf v ) = -3.9 
and —5.2 erg s" 1 cm ~ 2 at 1.2 Ryd. The dashed curves 
are the SEDs of T = 50, 000 (O-type) stars for flux lev- 
els \og(vf v ) = -2.4 and -4.4 erg s _1 cm ~ 2 at 1.2 Ryd. 
These normalizations were chosen because the lower flux 



grids illustrate minor modifications and the higher flux 
grids represent significant modifications to the UVB mod- 
els. The lower left panel shows the UVB photoionization 
grid superimposed upon the data. The mostly vertical 
curves are constant logarithmic neutral hydrogen column 
density, logA(Hi), and the mostly horizontal curves are 
constant logarithmic ionization parameter, log U. The col- 
umn density contours start at logA(Hi) = 20 cm~ 2 and 
decrease by 0.5 dex intervals toward the left. The ioniza- 
tion parameter contours start at logU = — 5 and increase 
by 0.5 dex intervals downward. The four remaining panels 
show the grids for the stellar ionizing continua. The grid 
contours follow the same pattern. 

The UVB model is consistent with a larger major- 
ity of the measured A(Fell)/A(Mgn) ratios. There are, 
however, several clouds that have measured ratios of 
log A(Fell)/A(Mgn) > and lie above the model grid. 
Adjusting the a-group to iron-group abundance pattern, 
[a/Fe], in the models is equivalent to "sliding" the model 
grid parallel to the A(Fen)/A(Mgn) axis. However, as- 
tronomical ratios range from [a/Fe] = (solar) for —1 < 
Z < to [a/Fe] = +0.5 (a enhancement) for Z < — 1 
(e.g., Lauroesch et al. 1996). The model grid would slide 
downward 0.5 dex for the a enhancement pattern. Since 
we assumed a solar pattern for the models the grid cannot 
be elevated further. Thus, it is not possible to model these 
these clouds by assuming a standard non-solar [a/Fe]. It 
is also not possible to simply assume a standard dust de- 
pletion pattern. In warm interstellar medium clouds, iron 
depletes more readily than does magnesium; iron and mag- 
nesium have depletion factors of S = —1.4 and —0.7, re- 
spectively (Lauroesch et al. 1996). Assuming a standard 
dust depletion would slide the model grid downward, ag- 
gravating the problem. 

Some of the elevated column density ratios could be 
an artifact of the VP fitting parameterization. At 
Keck/HIRES resolution, an unresolved Mgn line satu- 
rates (core intensity of zero) at log A ~ 13 cm~ 2 . In this 
regime of A(Mgn), the VP fitting simulations of blended 
absorption lines yielded roughly 10% by number VP com- 
ponents with similarly elevated A(Feii)/A(Mgn) ratios. 
It is thus plausible that a fair fraction of those clouds with 
log A ~ 13 cm~ 2 have elevated A(Feii)/A(Mgn) ratios 
as an artifact of the VP fitting process. 

However, this explanation cannot be invoked for those 
clouds with logA(Mgn) < 13 cm" 2 . The VP fitting 
simulations in this regime of A(Mgn) do not yield any 
components with elevated A(Feii)/A(Mgn) ratios. In- 
terestingly, in our sample, 15 clouds with logA(Mgn) < 
13 cm" 2 have A(Feii)/A(Mgn) > 0. It is worth not- 
ing that the ratios in these clouds hold some similarity to 
the "iron-rich" weak Mgn systems (Rigby et al. 2000). 
Though the global environments of the clouds in these 
strong systems may be very different than those of the 
weak systems, the physical processes that give rise to their 
abundance ratios (see discussion in Rigby et al. 2000) may 
be one and the same local to the clouds. 

Using Fen and Mgn as a diagnostic of the presence of 
stellar flux is not very useful because of the large degen- 
eracy between the column density ratios predicted for the 
different spectral shapes. The stellar SEDs tend to yield a 
maximum Fen to Mgll ratio that decreases as the stellar 
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flux is increased. This is especially true for the O stars. 
It is not ruled out that B stars can be contributing to the 
stellar flux for many of the model clouds illustrated. For 
B I stars, there would need to be ~ 40, 000 (100, 000) in 
number to make a minimum discernible contribution at a 
distance of R — 10 (20) kpc from the clouds (see Equa- 
tion 11 and Figure 11). 

5.4. Mgl to Mgll Ratio 

In Figure 13, we show results of the photoionization 
models for the Mgl to Mgll column density ratios over- 
plotted on the data. The photoionization grids are the 
same as described for the six panels of Figure 12. The 
Mgl balance is complicated by both a dielectronic recom- 
bination process and a charge exchange reaction with Hn. 
We compared model clouds for which these two mecha- 
nisms were included to those for which the mechanisms 
were "turned off" . In the relatively narrow density regime 
and temperature regime of the grids, the inclusion or ex- 
clusion of these processes did not influence the Mgl to 
Mgll ionization fraction. 

Illustrated in Figure 13 are two major discrepancies be- 
tween the data and the models. 

First, there is no single-phase photoionization 
model presented that can explain clouds having 
logiV(Mgi)/iV(Mgii) > -1.8. Though the results pre- 
sented here are statistical and involve comparing the data 
to a grid of models, we note that even detailed photoion- 
ization modeling involves struggles to reproduce observed 
Mgl to Mgll column density ratios (Churchill & Charlton 
1999; Ding et al. 2002; Rauch et al. 2002). 

Ding et al. (2002) argue that the majority of the Mgl 
gas arises in a lower temperature (T ~ 500 K), virtually 
non-turbulent, higher density phase than does the bulk 
of the Mgll absorption. That is, the Doppler b parame- 
ters of the Mgl clouds are, in reality, significantly smaller 
(b ~ 1 km s _1 ), than those of the bulk of the Mgll ab- 
sorption (b ~ 5 km s _1 ). If Ding et al. are correct, then 
"unphysical" observed large column density ratios are an 
artifact of VP fitting, which does not enforce any phys- 
ical connection between the Mgl and Mgll gas. In this 
scenario, much of the Mgll absorption would arise in a 
separate phase, warmer than that of the Mgl. If so, then 
the single-phase models presented in Figure 13 do not ap- 
ply to the a substantial fraction of the Mgl data. 

To avoid a high density solution, Rauch et al. (2002) 
suggest a large N(Mgi)/N(Mgu) ratio in the z — 0.5656 
system toward Q2237 + 0305 may indicate a large iV(Hl), 
on the order of a DLA. This cannot be the case for the 
majority of the large ratio clouds studied here because (1) 
it would imply too large a redshift path density for DLAs 
in this redshift range (Rao & Turnshck 2000), and (2) only 
three systems (S2, S5, and S7) are known to have near- 
DLA or DLA hydrogen column densities (Churchill et al. 
2000a,b). 

The second discrepancy between the data and the 
models in Figure 13 is that no model is con- 
sistent with the observed anti-correlation between 
log A f (Mgi)/7V(Mgn) and logiV(Mgii). This slope is 
Alog[iV(Mgi)/iV(Mgii)]/AlogiV(Mgii) = -0.73 over the 
range 11.5 < logTV(Mgn) < 14.5 cm~ 2 . 

We searched for other correlations and anti-correlation 



between the Mgl VP parameters and the Mgll VP 
parameters to see if the anti-correlation between 
log iV(Mgi)/iV(Mgii) and log AT (Mgll) could be induced 
by systematics in the VP fitting. If b(Mgl)/6(Mgll) > 1 for 
the majority of the clouds [an unphysical condition unless 
Mgl is associated with gas more turbulent than is Mgll], 
and/or if this ratio were found to increase as A^(Mgn) in- 
creases, this could induce the observed anti-correlation. 
However, accounting for uncertainties, all 6(Mgi)/6(Mgn) 
are consistent with unity and there is no correlation be- 
tween 6(Mgi)/6(Mgli) and A^Mgli). There is also no cor- 
relation between 6(Mgi)/6(Mgn) and A(Mgi)/A(Mgn). 
As such, there appears to be no obvious VP fitting correla- 
tion that is inducing the observed anti-correlation between 
logA r (Mgl)/A r (Mgn) and logJV(Mgii). 

Interestingly, the curves of constant ionization parame- 
ter for the highest log vf v stellar models run parallel to the 
slope, Alog[A^(Mgi)/A^(Mgii)]/AlogA r (Mgii), but they 
underprcdict the Mgl to Mgll column density ratio by 
roughly 1 dex. Also note, however, that the stellar fluxes 
for the illustrated models require an unreasonable number 
of ionizing stars; at a distance of 20 kpc, roughly 4 x 10 5 
O stars and roughly 2 x 10 6 B I stars are required. 

On the other hand, if the metallicity of the model clouds 
were increased about 1.5 dex, which would effectively move 
the photoionization grids directly to the right on the pan- 
els in Figure 13, then the model ratios could be consistent 
with the clouds having logTV(Mgn) > 13.5 cm~ 2 . How- 
ever, models of the lower Mgll column density clouds could 
not be made consistent with the data. We also note that 
the required metallicity would be 0.5 dex greater than so- 
lar (this argument is independent of abundance pattern). 

There appear to be two real discrepancies be- 
tween the observation and single-phase ionization mod- 
els: (1) single-phase models arc not consistent with 
logA r (Mgi)/A(Mgn) > -1.8; and (2) the models 
cannot explain the observed anti-correlation between 
log A r (Mgi)/A r (Mgn) and log iV (Mgll). As Ding et al. 
(2002) point out for two- phase models, each phase can 
have log A r (Mgi)/A^(Mgii) < -1.8 while the combined ab- 
sorption from these two phases can be made consistent 
with the total observed absorption strengths (equivalent 
widths) of Mgl and Mgll, respectively. Further explo- 
ration of the phenomenological and physical conditions 
governing the Mgl to Mgll absorption strengths is beyond 
the scope of this paper and is planned for future study. 

6. CONCLUSIONS 

We have performed an analysis on the Voigt profile fit- 
ting results reported by Churchill & Vogt (2001) for 23 
Mgn-selected quasar absorption line systems. The quasar 
spectra were observed with the HIRES instrument (Vogt 
et al. 1994) on the Keck I telescope and have a velocity 
resolution of ~ 6 km s _1 . For most of the Mgll systems, 
several Fen transitions and the Mgl A2852 transition were 
available for study. Further details on the data and the 
Voigt profile fitting have been presnted by Churchill (1997) 
and Churchill & Vogt (2001). 

In general, we have examined both the statistical prop- 
erties of the number of components, and the compo- 
nent column densities, Doppler parameters, and velocities. 
Where appropriate, we have compared our findings with 
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identical analyses on simulated data in order to ascertain 
possible statistical bias or systematics. These simulation 
results have provided insights into the degree to which 
an observational result is or is not an artifact of the Voigt 
profile fitting formalism. Further details of the simulations 
can be found in Churchill (1997). 

The statistical results of the VP analysis are summa- 
rized as follows: 

(1) The Fen and Mgn column densities are corre- 
lated at the 9 a level. A least-squares linear fit yielded 
log N (Fen) = 0.73 log JV(Mgll) + 3.0. A fit to the Mgi col- 
umn densities yielded log TV (Mgi) = 0.45 log iV(Mgll) + 
5.2. There is a 5 a anti-correlation of the Mgi column 
density with the ratio of the Mgi to Mgil column densi- 
ties. 

(2) The column density distributions for Mgil, Fell, 
and Mgi were parameterized with a power-law of the 
form f(N) oc N~ 5 . The maximum likelihood power-law 
indices were found to be S ~ 1.6, 1.7, and 2.0, respec- 
tively. Within 1 a uncertainties, the Fell and Mgil col- 
umn density indices are consistent with one another. How- 
ever, the Mgi index is significantly steeper than those for 
Fen and Mgil. There is a 5 a anti-correlation between 
the ratio A(Mgi)/A(Mgn) and A(Mgn) which follows 
logAT(Mgi)/iV(Mgii) = -0.73 log A(Mgn) + 7.6. 

(3) The modes of the Doppler parameter distributions 
were ~ 5 km s _1 for Mgil and Fen and ~ 7 km s _1 for 
Mgi. Based upon simulations, the true distribution modes 
could be 1-2 km s _1 smaller. The distribution tails at 
b > 10 km s _1 are usually due to blends in partially or 
fully saturated profiles. We found that the clouds are con- 
sistent with being thermally broadened, with most com- 
mon temperatures in the 30-40,000 K range. 

(4) We fitted a two-component Gaussian model to the 
two-point correlation function (TPCF) and found velocity 
dispersions of 54 km s _1 and 166 km s" 1 , with the nar- 
row component having twice the amplitude of the broader 
component. An analysis of the subsamplcs, demarcated 
by equivalent width, revealed that the clouds are dis- 
tributed more uniformly across the full velocity interval of 
an absorption profile as equivalent width increased above 
W r (2796) > 0.6 A. In the largest equivalent width sys- 
tems, extreme velocity clouds are less frequent than for 
systems with 0.3 < W r (2796) < 0.6 A. 

(5) A grid of photoionization models, in log[7 and 
log AT (Hi), were compared to the Fen to Mgil and Mgi 
to Mgil column density ratios. For the most part, 
the ratios are consistent with being photoionized by 
the ultraviolet extragalactic ionizing background. Stel- 
lar radiation has a difficult time producing clouds with 
log A r (Feil)/A r (Mgn) ~ 0. As the relative contribu- 
tion of the stellar flux increases, the maximum Fell to 
Mgil column density ratio allowed by the models de- 
creases. A small number of clouds are measured to have 
log 7V(Fell)/iV(Mgii) > 0, and there are no models that 
can produce this ratio. Similar conclusions hold for the 
Mgi to Mgil column density ratios. The majority of 
the clouds have log A r (Mgi)/A(Mgn) > -1.8, which can- 
not be made consistent with single-phase photoionization 
models. It is possible that the Mgi is arising in a sepa- 
rate, colder phase with a temperature of several hundred 
degrees. 



Here, we briefly take the opportunity to recap what 
has been observed and can be inferred about the absorber 
kinematics, geometry, and kinematic evolution. The Mgil 
absorbers break into kinematic subsystems, comprised of 
either a blend of components or a single component. In al- 
most all systems there is a single dominant subsystem with 
W / r (2796) ~ 0.2 A (on average). There are no examples 
of systems having two dominant subsystems with a sub- 
stantial velocity spearation. The dominant subsystem are 
comprised of clouds that are close in velocity and result in 
blended profiles. Kinematically, the "minor" subsystems 
either all lie to positive velocities or all lie toward negative 
velocities with respect to the dominant subsystem. These 
facts are suggestive of systematic kinematics and possibly 
a planar geometry for the dominant subsystem. 

For the largest equivalent width systems (sample E), 
there appear to be two distinct absorption profile mor- 
phologies (i.e. optical depth as a function of velocity). 
These systems have been classified as either "double" sys- 
tems or DLA/Hi-rich systems (Churchill et al. 2000b). 
The double systems are characterized by a large kinematics 
spread (greater than 200 km s _1 ) and a high variable opti- 
cal depth across the majority of the profile (double systems 
with smaller equivalent widths are characterized by weak 
absorption at the larger velocities). The DLA/Hi-rich 
systems are characterized by a moderate velocity spread 
(~ 60 km s^ 1 ) with and a black-bottomed, saturated pro- 
file. 

Now, also consider the following two facts: (1) The pop- 
ulation of the largest systems, with W,.(2796) > 1.0 A, 
evolves away the most rapidly compared to all Mgil sys- 
tems in the redshift range z = 2 to z = 0.3 (Stcidcl & 
Sargent 1992); and (2) The redshift number density of 
DLA systems is apparently not evolving over the studied 
redshift range (Rao & Turnshck 2000) . Since the DLAs ap- 
pear to have a distinct kinematics characterized by black- 
bottomed saturation and a moderate velocity spread, we 
can infer that this particular Mgil kinematic morphology 
is not evolving. Thus, it must be the double systems that 
are evolving away with redshift. We further note that 
it must be the subset of double systems, like SI and S9, 
which are characterized by large variations in optical depth 
across the majority of their velocity spread, that are dom- 
inating the evolution. This is because the weak higher 
velocity absorption in double systems, like S12, do not 
contribute significantly to the total equivalent width of 
the system. Evolution in these weak subsystems at higher 
velocities cannot be contributing to the observed evolu- 
tion. 

Thus, there appears to be an evolving population of very 
strong Mgil absorbers with complex kinematics that does 
not include DLAs. They must be Lyman-limit systems. It 
is an interesting question to ask if these systems are con- 
nected with the star formation history of their host galax- 
ies. Is there a population of post-bursting galaxies that 
fade at z <~ 2 but generate kinematically complex, highly 
spread out Mgil absorbers that passively evolve away to 
z ~ 1? Or, do these systems trace the "smooth" accretion 
of enriched gas from the intergalactic medium onto galax- 
ies at higher redshifts that becomes less dominant at lower 
redshifts (e.g., Murali et al. 2002). 
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Table 1 
Samples 



QSO 


Zabs 


ID 


Sample B 


Sample C 


Sample D 


Sample E 








0.3-0.6 A 


0.6 - 1.0 A 


> 0.6 A 


> 1.0 A 


0002 + 051 


0.85139 


SI 






X 


X 


0117 + 213 


0.57640 


S2 




X 


X 




0117 + 213 


1.04797 


S3 


X 








0420 - 014 


0.63300 


S4 




X 


X 




0454 + 039 


0.85957 


S5 






X 


X 


0454 + 039 


1.15325 


S6 


X 








0454 - 220 


0.47441 


S7 






X 


X 


0454 - 220 


0.48334 


S8 


X 








0823 - 223 


0.91102 


S9 






X 


X 


1101 - 264 


0.35900 


S10 


X 








1148 + 384 


0.55336 


Sll 




X 


X 




1206 + 459 


0.92760 


S12 




X 


X 




1222 + 228 


0.66805 


S13 


X 








1241 + 176 


0.55048 


S14 


X 








1248 + 401 


0.77296 


S15 




X 


X 




1254 + 044 


0.51939 


S16 


X 








1254 + 044 


0.93423 


S17 


X 








1317 + 274 


0.66005 


S18 


X 








1421 + 331 


0.90287 


S19 






X 


X 


1421 + 331 


1.17261 


S20 


X 








1634 + 706 


0.99024 


S21 


X 








2128- 123 


0.42973 


S22 


X 








2145 + 064 


0.79078 


S23 


X 









Table 2 



TPCF Parameters 



Sample 


Ax 






CT2 


A x 


ox 


A 2 /Ax 


(T2 




(xl0~ 3 ) 


km s -1 




km s _1 


(xl0~ 3 ) 


km s _1 




km s" 1 






With Doubles 






Without Doubles 




Full 


1.25 


54 


0.45 


166 


1.80 


60 


0.20 


151 


B 


1.51 


27 


0.83 


93 


1.51 


27 


0.83 


93 


C 


1.18 


55 


0.33 


273 


1.46 


58 


0.23 


195 


D 


0.91 


66 


0.55 


201 


1.66 


81 






E 


0.78 


67 


0.77 


154 


1.77 


84 
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Sample B W r (2796)<0.6 A 




Velocity km s 



Sample C 0.6^W r (2796)<1 .0 A 




-100 100 100 200 -200 -100 -500 -400 -300 -200 -100 100 200 300 -100 100 

Velocity km s 1 

Sample E W r (2796)^1 .0 A 




-100 100 -200 -100 100 200 300 400 500 600 -100 100 200 -100 100 -100 100 



Velocity km s 



Fig. 1. — Normalized HIRES/Kcck spectra of the Mgn A2796 absorption profiles presented in the system rest-frame velocity. The solid 
curves through the data (histogram) are the model spectra from Voigt profile (VP) fitting using the code MINFIT. The ticks above each 
continuum give the VP velocity centroids. The systems are presented by sample membership in order of increasing rest-frame equivalent 
width. Each system identification (see Table f ) and the equivalent width are given in the bottom part of the panels. 
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Fig. 2. — An illustration of the Voigt profile (VP) modeling of the Mgll AA2796, 2803 absorption profiles for the z = 0.7730 system in the 
spectrum of PG 1248 + 401 (ID S15). The solid curves through the data (histogram) are the model spectra from VP fitting using the code 
MINFIT. The ticks above each continuum give the VP velocity centroids. The individual VP components are shown as dotted curves. The 
"high velocity" cloud in this systems at v = +225 km s _1 is not shown. 




5 10 15 20 5 10 15 20 

N VP N VP 

Fig. 3. — (a) The distribution of VP components. The solid-line histogram represents the entire sample, whereas the shaded histogram 
represents the same sample but with the DLA and Hi— rich systems removed. The average number of clouds per absorber is 7.7 for both 
distributions. — (b) The rest-frame Mgll A2796 equivalent width versus the number of VP components. Two least squares fits are shown. 
The dotted line includes the full sample; the solid line is for the the same sample but with DLA and Hi-rich systems removed. See text for 
details. 
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log N(Mgll), cm" 
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log N(Mgll), cm" 



14 
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Fig. 4. — (a) The logarithmic VP column densities (in atoms cm -2 ) for Fell as a function of MgH. Upper limits for Fen are shown as 
open data points. The best— fit slope, shown as a solid line, is 0.73 ± 0.06. The reduced x 2 °f 30 provides a measure of the scatter about this 
relationship. The standard deviation of the fit is shown by the two dotted lines. — (b) The logarithmic VP column densities (in atoms cm -2 ) 
for Mgl as a function of Mgo. Upper limits for Mgl arc shown as open data points. The best— fit slope, shown as a solid line, is 0.45 ± 0.05. 
The reduced \ 2 °f 27 provides a measure of the scatter about this relationship. The standard deviation of the fit is shown by the two dotted 
lines. 
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Fig. 5. — (a) The logarithmic VP column density ratio of Fell to Mgn as a function of Mgn. The error bars are not shown for presentation 
purposes; errors and data point types can be viewed in Figure 4a. — (b) The logarithmic VP column density ratio of Mgl to Mgn as a 
function of Mgn. The best-fit slope, shown as a solid line, is —0.73 ±0.06. The standard deviation of the fit is shown by the two dotted lines. 
The error bars are not shown for presentation purposes; errors and data point types can be viewed in Figure 4b. 
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Fig. 6. — The logarithmic VP column density distribution functions for Mgn (top), Fen (center), and Mgl (bottom). The maximum 
likelihood results (unbinned data) for the relation f(N) oc N~ s arc 6 = 1.59 ± 0.05, 1.69 ± 0.07, and 2.02 ± 0.03 for Mgn, Fell, and Mgl, 
respectively. The number of clouds, m, in each sample is written in each panel. The shaded areas show the column density ranges where line 
blending leads to "partial completeness" (see text). 
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1 i i i r 




10 20 30 



b, km s 



Fig. 7. — The binned distributions of the VP Doppler b parameters for Mgn (top), Fell (center), and Mgl (bottom). The median Dopplcr 
parameters and the standard deviations of the distributions are (b) = 5.4 ± 4.3, 5.1 ± 4.1, and 7.7 ± 5.1, for Mgn, Fell, and Mgl, respectively. 
The number of clouds, m, in each sample is written in each panel. 



20 



CHURCHILL, VOGT, & CHARLTON. 




Fig. 8. — The thermal Doppler b component versus the totel Doppler b for Mgll. The deconvolution is based upon the Fen total Dopplcr 
b. The corresponding temperature (in Kelvin) is provided along the right hand axis. 
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Full Sample (W r >0.3A) 
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Fig. 9. — (top) The VP component velocity two-point correlation function (TPCF) for the full sample. The distribution is well described 
by a double Gaussian function with a narrow dispersion of <ri = 54 km s and a broader dispersion of <T2 = 166 km s^ 1 . — (bottom) The 
TPCF for the with double systems removed. The narrow dispersion is <ri = 60 km s~ 1 and the broader dispersion is <T2 = 151 km s _1 . 
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Fig. 10. — (a,b,c) The logarithmic VP column densities of Mgn, Fen, and Mgl as a function of cloud velocity. — (d,e) The logarithmic 
VP column density ratios of Fell and Mgl to Mgn as a function of cloud velocity. For all panels, the vertical scales are adjusted for ease of 
comparison. The error bars are not shown for presentation purposes; however, the errors and data point types can be viewed in Figure 4. 
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Fig. 11. — The minimum number of stars required at z = 1 to dominate the extragalactic ultraviolet background radiation as a source of 
photoionization. The curves are 0.4 dex lower at z = 0.5. Three types of stars are shown: O stars (solid), B I stars (short dash), and B V 
stars (long dash). 
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Fig. 12. — (upper left) Spectral energy distributions (SEDs) of the ionizing flux for the photoionization models. The thick solid curve is 
the Haard & Madau UV background for 2 = 1. The thin solid curves are T = 30,000 K Kuruzc stellar models (B stars) with normalization 
log(i//„) = —3.9 and —5.2 erg s _1 cm -2 . The thick dashed curves are T = 50,000 K Kuruzc stellar models (O stars) with normalization 
\og{vfi,) = —2.4 and —4.4 erg s _1 cm~ 2 . The photoionization model grids are shown in the remaining panels for log7V(Fell)/A r (Mgll) vs. 
log N (Mgn) as labeled with the data points overplotted. The error bars are not shown for presentation purposes; however, the errors and data 
point types can be viewed in Figure 4. The vertical curves are constant neutral hydrogen column density starting with logTV(Hl) = 20 cm -2 
on the right and decreasing by 0.5 dex moving leftward on the diagrams. The horizontal curves are constant ionization parameter, starting 
with log U = — 5 at the top and decreasing by 0.5 dex moving downward on the diagrams. (See text for details.) 




Fig. 13. — The same as for Figure 12 but for the logarithmic VP column density ratio of Mgl to Mgo as a function of Mgll. (See text for 
details.) 



